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1 Introduction 

The traditional algorithm of collisional process in DSMC method, used for 
calculation of velocities of molecules after collision, is based on the conser- 
vation laws for linear momentum and total energy . This algorithm does 
not provide conservation of angular momentum with respect to some axis, 
but for non-rotating flow this circumstance is insignificant, since the mean 
value of angular momentum of such kind of flow is equal to zero. 
The situation changes for axisymmetric flow with rotation. In this case 
the traditional procedure leads to some internal source or sink of angular 
momentum that may distort the flowfield. 

In the paper the discussed effect is recognized and investigated. The angular 
momentum conservative (AMC) algorithm of collisional process is proposed. 
The draft version of the paper is published in |2] . 

2 Problem Formulation 

The investigations were made for one-dimensional axisymmetric flow of 
monoatomic gas in the tube with specular wall. The evolution of initially 
swirling flow during the time was studied. The gas in the tube with uniform 
initial density no and temperature To and equilibrium distribution function 
was assumed to start rotating as a solid body, the initial tangential velocity 
on the wall being equal to v v . For given molecular model the problem con- 
tains two governing parameters: the Knudsen number Kn, defined by the 
ratio of the mean free path lo in the gas at t = to the radius of the tube 
r t 

Kn = — 

n 



and the speed ratio 



Sip = v v /y / 2RT 



The simulations were made on uniform grid with 10 4 simulated molecules. 
Standard NTC DSMC procedure [T] was applied. 

The initial values of total mass Mq, angular momentum Iq and energy Eq 
of the gas may be defined by the relations 

Mq = m Nq — mriQiTrl (1) 



N 



Ia = 171 2^ Vi r * ( 2 ) 

i=l 
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E ° - m E -f ( 3 ) 

i=l 

where Nq is the total number of simulated molecules, m - the mass of 
molecule, V - the velocity vector, v - the tangential component of this 
vector. By the same relations may be determined the total energy E and 
total angular momentum / of the gas during the time. 

For any values of Kn and S v there should be no temporal dependence of E 
and /. The computations show that this condition is satisfied only for E. 
FigGJilluatrates the dependencies of I/Iq on the time for some of considered 
variants for S v = 1. The time is normalized by the value of to — 2 7rr t /i> v , 
so t/to represents the number of revolutions of solid body with v(r t ) = v v . 
As can be seen from Fig^ for all of the variants the decreasing of I/Iq with 
increasing of t/to is observed. For given Kn the effect is more pronounced 
for solid sphere molecules (SS) in comparison with VSS molecular model 
for Maxwell molecules (MM). The effect depends on the Knudsen number: 
the less Kn the more the effect. But the highest sensitivity of the effect 
is observed to the grid size: changing the number N of cells from 50 to 25 
causes the same effect as the decrease of the Knudsen number approximately 
by 4 times. 

Due to the dependence of I(t) it is difficult to obtain the steady solution 
of the problem. The dissipation of rotational contribution of energy into 
the heat takes place during the time. The only steady solution that may 
be obtained with traditional collisional procedure is terminal one, when at 
t — > oo the gas will stop rotating, having uniform density no and getting 
warm to the temperature 

T = 2Eq/3RM 



3 AMC Algorithm 

The reason of temporal dependence of total angular momentum is quite 
clear: the colliding molecules have different radial coordinates r\,r^, that 
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Figure 1: Temporal dependence of total angular momentum 



leads to the difference in angular momentum 

ni\ Vi ri + nii "2 r 2 

of these molecules before and after collision (v denotes the tangential com- 
ponent of velocity of the molecule) . To provide the conservation of angular 
momentum the algorithm of collisional process should be changed. One 
possible variant of AMC algorithm is as follows. 

Let us denote by u the vector containing axial and radial components of the 
velocity of molecule. The post-collision velocities calculated by traditional 
procedure will be denoted by symbol *, while the symbol ' will be used for 
these velocities in the proposed algorithm. 

The algorithm is based on conservation law of angular momentum mvr for 
tangential component of velocity instead of linear momentum mv. That is 
why the post-collision tangential velocities are calculated by usual relations 
PQ, with the value of m, being replaced by the product m, n- The obtained 
velocities v[ , v' 2 provide precise conservation of angular momentum, but vio- 
late the energy conservation. To provide energy conservation the correction 
of velocities u\ and u\ is needed. To do this correction the analysis of 
energy defects 

dEi = nii (vf - v* 2 ) 
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should be made. Three possible relations between energy defects dEi and 
velocities u* are as follows. 



1. If 

and 



mi u* 2 > dE\ 



m 2 u 2 2 > dE 2 

both components of vectors u* should be corrected by the factors Ci 



= \ I - dEi/rriiU 



*2 



2. If 

mi u* 2 + m 2 u*2 > dE\ + dEi 

the correction factor c for components of both velocities u* should be the 
same and equal to 



1 - {dEx + dE 2 )/(m lU f + m 2 u* 2 ) 



3. If none of the above conditions is satisfied, the collision is considered to be 
"bad" and is not performed, i. e. the molecules conserve their pre-collision 
velocities. 

The computations show that the relative number of collisions of types 1 
and 2 is about 99% and 1%, respectively. The relative number of "bad" 
collisions never exceeds 10 -4 , so the effect caused by the neglecting of these 
collisions is small. 



4 Steady Solution 

The described algorithm provides precise conservation of angular momen- 
tum and energy. All the curves, presented in Fig ^transform into straight 
line I /Iq = 1, if this algorithm is applied. 

AMC algorithm enables to obtain the steady solution of the considered 
problem. This solution proves to be the same, as the prediction of Navier- 
Stokes theory and has the following features. 

1. The flow is isothermal 

T = const (4) 

2. The gas rotates as a solid body 

v = v w r/r t (5) 
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3. The radial distribution of density is described by the relation 



n(r) = n(0) exp(u 2 /2i?T) (6) 

4. The solution is completely determined by the values of Mq, Iq, Eq and 
does not depend on the way of initial swirling. 

5. The solution does not depend either on the molecular model or on the 
Knudsen number. 

These features of the considered flow are important for clarifying the nature 
of the Ranque effect 0] . 

Three unknown values v w ,T, n(0) determining the radial distribution of 
parameters of steady solution may be defined based on the values of Mg, Iq 
and Eq from the relations 



7rr>n(0) (eXp g 1} = M (7) 



MoV - Tt ^(exp^-1) = 7 ° (8) 

^TMo + ^»=B„ (9) 
2 2r t 

where speed ratio S w of the flow at r = r t is defined by 



S w = v w /V2RT (10) 

The left hand sides of the relations 0-03 representing the total mass, angular 
momentum and energy of the flow may be obtained based on the relations 
El-El by simple integration. 

The comparison between numerical and analytical radial distribution of 
parameters of steady flow for VSS molecules, Kn — 0.1, 5^ = 1, N = 50 
is made on Figs. [5J-0] for density, tangential velocity and temperature, 
respectively. The numerical results are shown by solid circles while solid 
lines represent the results of analytical solution [3-13 As can be seen from 
these Figures, numerical and analytical results are in excellent agreement. 
For Kn = 10 the steady solution was found to be the same as for Kn = 0.1. 

It should be noted that the steady flow in the considered problem is charac- 
terized by locally Maxwellian distribution function for any molecular models 
and Kn < oo. 

AMC algorithm of collisional process may be applied for simulation of non- 
rotating flows too. The experience of such application to the flow in gas- 
dynamic window |S| shows that the results obtained with proposed and 
traditional algorithms are the same within the expected statistical scatter. 
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Figure 2: Radial distribution of density 
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Figure 3: Radial distribution of tangential velocity 



To stay the effect of non-conservative nature of traditional collisional proce- 
dure in the real conditions, some test computations of considered flow were 
made by AMC algorithm for adiabatic tube wall with diffuse-specular scat- 
tering of molecules for different values of accommodation coefficient a. FigEl 
illustrates the dependencies of I/Io on the time for VSS molecules, Kn = 
0.1, = 1, N = 50 and 4 values of a in the range a = 2 ■ 10~ 3 - 10~ 4 . 
The comparison of these dependencies with those presented in Fig^ reveals 
that the traditional collisional procedure causes the same change in angular 
momentum as the diffuse-specular wall with the accommodation coefficient 
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Figure 4: Radial distribution of temperature 




Figure 5: Temporal dependence of total angular momentum 
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5 Conclusion 



The proposed algorithm of collisional process provides precise conservation 
of angular momentum and energy and may be used for direct simulation of 
flows with rotation. 

It should be noted, however, that the improvement of the results, obtained 
by this algorithm may be important only for the problems without external 
source or sink of angular momentum. If solid surface with diffuse reflection 
is present in the flowfield, the changes in angular momentum, caused by this 
surface will be much greater than those, caused by non-conservative nature 
of traditional algorithm of collisional process. For such kind of flows there 
is no need to employ AMC algorithm, especially taking into account that it 
takes 10% - 15% more CPU time. 

Nevertheless, for some problems the application of AMC algorithm is advis- 
able. An example of such kind of problem is the Ranque effect which 
nature is determined by conservation of angular momentum of the gas flow 
in the vortex tube 0] . 

References 

[1] Bird G.A., Molecular Gas Dynamics and Direct Simulation of Gas 
Flows, Oxford, Clarendon, 1994. 

[2] Skovorodko P. A., Angular Momentum Conservative Algorithm of Col- 
lisional Process in DSMC Method, Rarefied Gas Dynamics, Book of 
Abstracts of 21st Intern. Symp., Marseille, France, 1998, Vol. 2, pp. 
107 - 108. 

[3] Ranque G.J., Experiences sur la Detente Giratoire avec Productions 
Simultanees d'un Enchappement d'Air Chaud et d'un Enchappement 
d'Air Froid, J. Phys. Radium, 4, pp. 112 - 114, 1933. 

[4] Rebrov A.K. and Skovorodko P. A., The Ranque Effect in Rarefied Gas, 
Rarefied Gas Dynamics, Proc. 21st Intern. Symp., Edited by R. Brun, 
R. Campargue, R. Gatignol and J.-C. Lengrand, Cepadues-Editions, 
Toulouse, France, 1999, Vol. II, pp. 221 - 228. 

[5] Karsten V.M., Sharafutdinov R.G., Skovorodko P.A. and Skrynnikov 
A. V., Experimental Investigation and Direct Simulation of Flow in 
Gasdynamic Window, ibid., Vol. I, pp. 655 - 662. 



8 



